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ABSTRACT- When a flat stress- free surface (i.e., the ground in seismological applications) sep- 
arating air from a isotropic, homogeneous or horizontally-layered, solid substratum is solicited 
by a SH plane body wave incident in the substratum, the response in the substratum is a sin- 
gle specularly-reflected body wave. When the stress-free condition, equivalent to vanishing surface 
impedance, is relaxed by the introduction of a spatially- constant, non- vanishing surface impedance, 

■ the response in the substratum is again a single reflected body wave whose amplitude is less than 
the one in the situation of a stress-free ground. When the stress- free condition is relaxed by the 

i-H | introduction of a a spatially-modulated surface impedance, which simulates the action of an un- 

even (i.e., not entirely- flat) ground , the frequency-domain response takes the form of a spectrum 
q . of plane body waves and surface waves and resonances are produced at the frequencies of which one 

*5a ! or several surface wave amplitudes can become large. It is shown, that at resonance, the amplitude 

of one, or of several, components of the motion on the surface can be amplified with respect to the 
situation in which the surface impedance is either constant or vanishes. Also, when the solicitation 
'— 'I is pulse-like, the integrated time history of the square of surface displacement and of the square of 

' velocity can be larger, and the duration of the signal can be considerably longer, for a spatially- 

■ modulated impedance surface than for a constant, or vanishing, impedance surface. 
Os ■ 

co ; 

(N . 

q '. 1 Introduction 

: 

An important question in seismology, civil engineering, urban planning, and natural disaster risk 
& ', assessment is: to what extent does surface topography of different length and height scales (ranging 

. — ■ from those of mountains and hills to city blocks and buildings) modify the seismic response (in terms 

^ of cumulative motion intensity and duration) on the ground? 

^ There exists experimental evidence (Singh and Ordaz, 1993; Davis and West, 1973; Griffiths 

and Bollinger, 1979) that this modification is real and can attain considerable proportions. 

Some theoretical studies (Wirgin, 1989; Wirgin 1990; Wirgin and Kouoh-Bille, 1993; Groby, 

^ . 2005) seem to indicate that such effects are indeed possible, but various numerical studies ( Bouchon, 

1973; Bard, 1982; Sanchez-Sesma, 1987, Geli et al., 1988; Wirgin and Bard, 1996; Gueguen, 2000; 
Clouteau and Aubry, 2001; Gueguen et al., 2002; Semblat et al., 2003; Tsogka and Wirgin, 2003; 
Boutin and Roussillon, 2004; Kham, 2004; Groby and Tsogka, 2005) yield conflicting results in that 
some of these point to amplification, while others to very weak effects, or even to de-amplification. 
Contradictory results are also obtained regarding the duration of the earthquakes. 



1.1 Sites 



In figs. Hand El we give examples of natural and man-made sites respectively with quasi-periodic 
features that can be studied by the methods of this investigation. 




Figure 1: Appalachian mountains. 




Figure 2: Istanbul, Turkey. 
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Figure 3: Periodically- uneven ground. 

r o hH 




Figure 4: Periodically- modulated impedance flat ground. 
1.2 Physical configurations 

The (assumed) periodic uneveness of the ground (an example of which is given in fig. |3J) is accounted 
for by a suitably-chosen spatially-periodic flat surface impedance function (see fig. HJ. 



1.3 On the notion of impedance 

We shall employ the notion of (surface) impedance employed in the civil engineering (Gueguen, 
2000), (Roussillon, 2006) community wherein the three components of impedance are M, K, C 
which designate mass, stiffness and damping ( real constants) respectively. 

Thus, the mechanical impedance is complex and designated by Z = R — iX, where R = C is 
the "resistive" part, and X = uM — K/u the "reactive" part. 

The reactance is "inductive" if lo 2 > K/M (i.e., X > 0) and is "capacitive" if K/M (i.e., 
X < 0. 
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2 Governing equations 

2.1 Mathematical translation of the boundary value problem in the space- 
frequency domain 

n(x, u) + it,22 (x, w)] + puj 2 u(x,uj) = ; Vx € , (1) 
zcjZ(xi, w)it(x, cj) + /Uit, 2 (x, cj) = ; Vx € T , (2) 
-u d (x, u;) := u°(x,u;) — u J (x, oS) ~ outgoing waves ; ||x|| — > oo , x € , (3) 



uVx, w) = ^(w) exp[i(^xi - /c|:c 2 )] ; Vx G , k\ = k sin 6 l , k l 2 = kcos6\ k = —, (4) 

c 

wherein: u is the total displacement field, the (unknown) diffracted field, u l the (known) incident 
field , 9 l the angle of incidence with respect to the x 2 axis, A 1 (to) the incident pulse spectrum. 

Z(x\, uj) := R(xi, u) — iX(x\, u) , (5) 

When the impedance vanishes for all x±, the boundary condition becomes that of a flat, stress-free 
surface. 

Otherwise, the impedance boundary condition is supposed to simulate the presence of a topo- 
graphically- uneven stress-free surface. 

The periodic nature of Z is expressed by: 

Z(x\ + d,u) = Z(x\,lo) ; \/x\ £ R , (6) 



Z\ = J ^ Z(a;i) exp ; VZ G Z . (7) 

It is also assumed that i?(xi,w) > 0, i.e., the impedance is passive and dissipative. 

2.2 Diffracted field representation 

00 

i (x 1 ,x 2 )= ^2 ^nexp[i(/ci n xi + k 2n x 2 )] ; V(xi,x 2 ) G ft , (8) 

n=— 00 

2n7r 



fcin = k\ + — , k 2n = y Zc 2 - /c 2 n , Kfc 2n > ; 9/c 2 n > , cj > . (9) 

The diffracted field is a discrete sum of plane waves: those for which k 2n is real are propagative 
(or homogeneous) and those for which k 2n is imaginary are evanescent or (inhomogeneous). 
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2.3 Result of the introduction of the field representation into the boundary 
condition 

oo 

+ ln5nm]B n = -2Z m A l ; Vm € Z , (10) 

71= — OO 

wherein 

B n = A n -A i S n0 , in-=— ; VneZ. (11) 
The linear system can be written as the infinite-order matrix equation 

Ef = g. (12) 



3 On the possibility of anomalous fields in the general case of a 
periodic, passive, spatially non-constant, surface impedance 

A mode of the configuration is obtained by turning off the solicitation in the matrix equation, i.e., 
g = 0, wherein is the null vector. 

We are thus faced with the equation Ef = 0, whose solution is trivial (i.e., f = 0), unless 

det(E) = . (13) 

An " eigenvalue" is a value of k\ n for which the determinant vanishes at a given frequency. Another 
way of putting things is to fix k\ n and look for the frequencies ("natural frequencies") that lead to 
a vanishing determinant. 

When the configuration is such that one of the k± n is an eigenvalue, and the frequency is a 
natural frequency, then the system is said to be in a state of resonance. 

When this happens, the determinant of E is either small or nil, which means that the inverse 
of E is either large or infinite and that consequently one or more entries in the vector f of the 
scattered plane wave coefficients are either large or infinite. 

Consequently, we can expect (since the field is a discrete sum of scattered plane waves) that the 
field may become large at resonance (in the presence of not too much dissipation and/or radiation 
damping) . 

From now on, we call det(E) = the general dispersion relation. In the case of spatially- 
constant surface impedance Zq = only the zeroth-order diffracted wave is generated, so that the 
dispersion relation is 

C + ^ = 0. (14) 

U) 

and the latter has no solution for a passive impedance (i.e., $l( > 0) due to the fact that k\ is real. 

This means that no resonance can be produced for a constant, passive impedance surface. When 
the surface impedance function Z is not spatially-constant, it is much more difficult to obtain a 
meaningful expression of the dispersion relation. 

Some insight may be gained by turning to an iteration method of solution 



2AZ m ^ N ^ m Z m „ n B n m = - N ,.. A ..,N; 1 = 1,2, (15) 
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which suggests that B m can become large for 



Z + {^ = 0. (16) 

UJ 

This cannot occur for m = for the previously-mentioned reason. It can occur, if at all, only for 
m ^ 0. Recall that it was assumed that 3^=0, Qk = and 

Z(xi) = R{xi) - iX(xi) => Z m (x\) = Rm(xi) -iX m [x\) . (17) 

Consider the case of vanishing resistance, i.e., R{x\) = in which the approximate dispersion 
relation is 

-iX + ^k 2m = . (18) 
uj 

Since uj > 0, > 0, ffltzm > 0, and 5?/c2m > 0, the second term in the previous equation is either 
positive real (for real kirn) or positive imaginary (for imaginary k2 m ), so that the sum of the two 
terms can vanish only if 

X > and m 2m = . (19) 

The first of these requirements means that the impedance must be inductive for it to be possible to 
obtain resonant behavior. Thus, a possible explanation of why several researchers, who employed 
the impedance concept to account for ground uneveness, have not been able to obtain anomalous 
fields is that their impedance functions were such that Xq < in the frequency range of the 
incident pulse. The second requirement, i.e., 3?&2m = 0, means that resonances can occur only for 
the evanescent waves in the plane wave representation of the scattered field. Thus, we can expect 
the amplitude of the m-th order evanescent wave to become infinite (for R = 0) or large (for R > 0) 
at resonance, which is another way of saying that a surface wave (evanescent waves are of this sort) 
is strongly excited at resonance (all the more so the smaller is R). 

This picture is only partially true, because the approximate dispersion relation may account 
only poorly for all the features of the solutions of the general dispersion relation, which is the case 
if the uneveness of the ground is large. 

The fact that B m becomes large at resonance is a necessary condition for the ground motion 
to be large at this resonance frequency, but is not a sufficient condition due to the fact that the 
diffracted field is composed not of one evanescent plane wave, but of a sum of both propagative 
and evanescent plane waves, and this sum can be of modest proportions even when one of its 
components is large. Such modest fields at resonance are due mainly to radiation damping. 

3.1 The case of spatially- sinusoidal surface impedance 

For the computations, wee shall choose: 



Z(xi,u) = CM 



. 27r 

1 + h COS —rX\ 

a 



(20) 



When h = 0, we re-encounter the case of a constant impedance on flat ground, so that h is a 
measure of the uneveness of the ground. 
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4 Computations 



We shall be interested in the following quantities indicative of possible anomalous effects provoked 
by the uneveness of the ground: 

• the duration of significative seismic ground motion, 

• the peak value of the ground displacement, 



the amplification factors (i.e., amplification if the factor > 1, de-amplification if the factor < 1) 
of the time integral (over [0, r]) of the squared displacement u (velocity v) for the modulated 
impedance surface relative to the time integral (over [0, r]) of squared displacement (velocity) 
for the h = impedance surface: 

_ j-[u(0At\C,h)?dt _ £[v{o,o,t\c,h)] 2 dt 

x[( » n > T) £[u(o,o,t\c,o)] 2 dt ' v[( " n ' T) / >(o,o,t|c,o)pdi ' [ZL) 

the amplification factors of the time integral (over [0, r])of squared displacement (velocity) 
for the modulated impedance surface relative to the time integral (over [0, r]) of squared 
displacement (velocity) for the stress-free surface: 

_ £K0 1 MC 1 ^ _ / >(o,o,*|C,/Q] 2 <ft 

m ' ,T) ~ Qu(0,0,t\0M 2 dt ' U ' ' J " / T [v(0,0,i|0,0)] 2 dt • 1 ' 

The amplitude spectrum of the incident plane wave is that of a Ricker pulse, i.e., A l (u) = 

The other parameters are: a = 1, (3 = 4, s l = 0.3, r = 0.1, m = 2, k = 1, s = 2. 

The results of the computations, peformed by direct resolution of a suitably-chosen finite dimen- 
sion version of the matrix equation, are given in figs. EJandlHlfor h = 0.1 and h = 1 respectively. 
The first line in these figures contains the spectra, the second line the time histories, and the 
remaining lines the ratios of cumulative response. 



5 Comments 

The spectra reveal the existence of a series of resonances, as predicted by the theory. 

The peak displacement (in the time domain) on a spatially-modulated impedance ground is 
always smaller than the peak displacement on the flat, stress-free ground. 

The maximum values of the ratios r/, Xi v an d v first increase and then decrease with increasing 
h, with the maximum amplifications (j) ~ 1.2, x ~ 11-7, v ~ 0.9, v ~ 15.7) being attained for 
h = 1.5. 

It is possible to obtain substantial amplification of the ratios with respect to the spatially- 
constant impedance flat surface and much less amplification or even deamplification of the ratios 
with respect to the stress-free flat surface. 

This shows that although there is no doubt that spatially-modulated impedance surfaces can 
give rise to beating and very long durations, attaining in some of these examples of the order of 2 
min versus of the order of 10-20 sec for zero or non-zero constant impedances), there is no clear-cut 
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Figure 5: h = 0.1 . Small-scale anomalous effects. 
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Figure 6: h = 1.5 . Large-scale anomalous effects 



answer to the question of whether spatial impedance modulations, which simulate the existence 
of ground uneveness, systematically result in amplification or even deamplification of cumulative 
ground motion displacement and velocity, especially with respect to motion of a stress-free flat 
ground. 
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